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Abstract 

We study the hadro- and photo-production of cc and bb mesons at 
low transverse momentum to high orders in the relative velocity of the 
pair, v, in non-relativistic QCD. We evaluate cross sections to order v 7 
for rjc, sufficient for studies of photo-production in the almost-elastic 
region. For all other charmonium states we find the cross section to 
order v 9 , sufficient for studies of the ratio of Xci and \c2 production 
rates. We find recurrence formulas for generating terms at even higher 
orders, should they be needed. 
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1 Introduction 



Recent progress in the understanding of cross sections for production of 
heavy quarkonium resonances has come through the non-relativistic QCD 
(NRQCD) reformulation of this problem [jlj. Although the production of 
the pair is dominated by short distance scales of order 1/m (where m is the 
heavy quark mass), longer non-perturbative scales play an important role. 
The projection on to a specific quarkonium state involves length scales such 
as 1/mv, l/mv 2 , etc., where v is the (dimensionless) velocity of either of the 
heavy fermions in the rest frame of the pair. 

Factorisation of the physics at the long and short scales has been proven 
in the NRQCD formalism for processes dominated by a large transverse mo- 
mentum [[|. The resulting cross sections are a double power series in the 
QCD coupling a s evaluated at the NRQCD factorisation scale /io, and the 
velocity v. Often, higher orders in v involve the previously neglected colour- 
octet states of the heavy quark pairs. For bottomonium states, v 2 <C a s {m 2 ) 
and hence colour octet contributions are often not very significant and the 
expansion is close to the normal perturbative expansion. For charmonium 
states, a numerical coincidence, v 2 ~ a s (m 2 ), makes the double expansion 
more complicated. 

The formalism has been successfully applied to large transverse momen- 
tum processes |§. Interestingly, inclusive production cross sections for char- 
monium at low energies, dominated by low transverse momenta, also seem to 
have a good phenomenological description in terms of this approach 0, y, g]. 
It was argued that a better understanding of such cross sections can be 
obtained if the higher order terms in v and a s are used. The argument is 
simple. Total inclusive J /if) cross sections arise either from direct J/ if) pro- 
duction (which starts at order a s v 7 ) or through decays of x states. Now Xo 
and X2 are first produced at order a s v 5 , whereas Xii which has the largest 
branching fraction into J /if), is produced only at order a s v 9 . Hence, a full 
understanding of these cross sections requires the NRQCD expansion upto 
order a s v 9 . We calculate these terms here. Of course, higher orders in a s 
at lower orders in v may be equally important. First attempts at computing 
these have been made 0j. 

The plan of this paper is the following. In section 2 we set out the nota- 
tion with a brief review of the threshold expansion technique ||, and some 
extensions, which we shall use for our computation. In the next section we 
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perform a Taylor expansion of the transition matrix elements in the relative 
momentum of the heavy quark pair and give recurrence relations for the Tay- 
lor coefficients to all orders. The cross sections are listed in section 4. We 
conclude in section 5 with a discussion of the phenomenological applications 
of this computation. The appendix contains a discussion of the complete 
specification of states and operators required for computations beyond the 
leading orders of NRQCD. 



2 The Threshold Expansion 

The NRQCD factorisation formula for inclusive production of heavy quarko- 
nium resonances H with 4-momentum P is 

1 P 

where $ is a flux factor. The coefficient function CV, is computable in per- 
turbative QCD and hence has an expansion in the strong coupling a s (evalu- 
ated at the NRQCD cutoff), whereas the matrix element is non-perturbative. 
However, in NRQCD, it has a fixed scaling dimension in the quark velocity 
v. Consequently, the cross section is a double power series in a s and v. 

The fermion bilinear operators /Q are built out of heavy quark fields 
sandwiching colour and spin matrices and the covariant derivative D. The 
specification of the composite labels % and j is given in appendix |A]. They 
include the colour index a, the spin quantum number S, the orbital angular 
momentum L (found by coupling the derivatives), the total angular momen- 
tum J and the helicity J z . At low orders in v this set is sufficient to fix the 
operators. Since the hadronic projection operator 

11(F) = -£\ H ^)(H,s\, (2.2) 

s 

(where s denotes hadronic states with energy less than the NRQCD cutoff), 
is diagonal in these quantum numbers, it is clear that the operators and 
JCj in eq. (|2.1|) are restricted to have equal L, S, J and J z . For a more 
detailed discussion see the appendix. 
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The J^-dependence of these matrix elements can be factored out using 
the Wigner-Eckart theorem — 

(icMml) = ^r[ "C s+1 L^ (2.3) 

where the first factor on the right comes from a Clebsch-Gordan coefficient. 
This factor is conventionally included in the coefficient function. In the 
reduced matrix element O h , we have introduced a new label iV which is 
the number of derivatives in each fermion bilinear (see the appendix). In 
agreement with the notation of we write for the off-diagonal operators 

(W(H)!C]) = ^ rT ^( 25+1 ^, 2S+1 Lf). (2.4) 

The power counting rule for the matrix elements in eq. Q2.1D is — 

d = 3 + N + N' + 2(E d + 2M d ), (2.5) 

where E d and M d are the number of colour electric and magnetic transitions 
required to connect the hadronic state to the state /Q|0). Note that at low 
orders in D, N = L, and the more familiar rules are obtained. An example 
is provided by the off-diagonal matrix element 

VT (X ^o 2 ) <Vxn(77 c ) X t ( _ * D) . (-lD)x) + h.c. (2.6) 

which scales as v 5 . The — l/y/3 factor on the right is a trivial Clebsch-Gordan 
factor, explained later. 

We choose to construct the coefficient functions using the "threshold ex- 
pansion" technique of ||. This consists of calculating, in perturbative QCD, 
the matrix element M. connecting the initial states to final states with a 
heavy quark- ant iquark pair (QQ), and Taylor expanding the result in the 
relative momentum of the pair, q, after performing a non-relativistic reduc- 
tion of the Dirac spinors. The resulting expression is squared and matched to 
the NRQCD formula of eq. (|2.1|) by inserting a perturbative projector onto 
a non-relativistic QQ state between the two spinor bilinears. The coefficient 
of this matrix element is the required coefficient function. 

Symbolically — 

J2\M\ 2 = £a,(/CTI(QQ)/Cj), (2.7) 

pol ij 
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where the left hand side is Taylor expanded in q. Each factor of q Fourier 
transforms into a factor of the covariant derivative D on the right hand side. 
Since each matrix element on the right of eq. Q2.7| ) corresponds to an unique 
matrix element in eq. ( [2.1|) , the order upto which the Taylor expansion is 
to be performed is determined by the scaling of the non-perturbative matrix 
elements with v. Since we require a classification of operators by the angular 
momentum, it turns out to be very convenient to use spherical tensor meth- 
ods. These were used in an earlier paper || and are used more extensively 
here. 

In this paper we evaluate the cross sections to order a s v 9 . The Taylor 
expansion order, N + N' < 6 is obtained by setting d = 9 and Ed = Md = 
in eq. ( |2.5| ). Furthermore, since we examine the leading term in perturbation 
theory, the perturbative projector has only one term — 

n(QQ) = \QQ){QQ\. (2.8) 

In agreement with M we use the relativistic normalisation 

(Q(p,OQ(q,V)\Q(p',e)Q(q , ,v')) = ±E p E q (27r) 6 8 3 ( P -p')8 3 (q-q'), (2.9) 

with the spinor normalisations £t£ = rfrj = 1. Since E p = E q = \Jm 2 + q' 2 , 
expanding this in q 2 allows us to write the spinor bilinears in terms of tran- 
sition operators built out of the heavy quark field. For example, 

?v = ^-(QQ(q)\^x\o) - ^(QQ(q)\^v ■ d x |o) + • • • (2.10) 

Conventionally the coefficient functions and matrix elements in eq. ( |2 . 1| ) 
were written with a non-relativistic normalisation of the hadron states. In 
the threshold expansion technique it is more convenient to retain a relativistic 
normalisation similiar to that in eq. ( |2.9| ). The result for the cross section is 
the same in either case, since a change in the definition of the matrix element 
is compensated by a change in the coefficient function. To leading order in 
q, the matrix elements in the notation of || have to be multiplied by 4m 
to obtain those in the relativistic normalisation ||. As higher orders the 
relation is more complicated. In this paper, we shall work entirely with the 
latter. 



4 



3 The Matrix Elements 



To leading order in a s , the kinematics is very simple. The momenta of the 
initial particles are p\ and p 2 - We take pi to lie in the positive ^-direction 
and pi to be oppositely directed. The momentum of the meson, P = pi +P2- 
As a result, s = P 2 = M 2 , where M is the meson mass. 

The 4-momenta of Q and Q (p and p respectively) are written as 

p = ^P + Lrf and p = l -P-LflK (3.1) 

Note that p 2 = p 2 = m 2 , where m is the mass of the heavy quark. The 
space-like vector q is always defined in the rest frame of the pair, and Lj 
boosts it to any frame. We shall use Greek indices for Lorentz tensors and 
Latin indices for Euclidean 3-tensors. 

The following relations are easy to prove — 

M , M 

Px-Lj = — -z h and p 2 ■ Lj = —Zj, (3.2) 

where z is the unit 3-vector in the z-direction. They are consistent with the 
identity P ■ Lj — 0. We shall use the two identities || 

Lj ■ L k = -5 ij7 and Me ijk L ak = e^L^P? . (3.3) 

Other relations can be written down 0, but are not important for our com- 
putations. Note our convention e 0123 = 1. 

This technique also depends on the usual non-relativistic reduction of 
Dirac spinors which gives rise to the identities 



u(p)n 5 v(P) = TjHrP^V - 2iL^e mnj q n ^^r]. 



(3.4) 



Here £ and rj are Pauli spinors and <r J are the usual Pauli matrices. M 
(M 2 = 4m 2 + 4g 2 ) is the invariant mass of the QQ system. 

We work in a class of ghost-free gauges called the planar gauges [|10| 
These are defined by the polarisation sum for gluons 



E C »< A (P) = d M = -9»v + —y(P t yv+PvV tl ). (3.5) 
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The propagator for a gluon of momentum p is then given by G^ v (p) = 
dfj,u(p)/p 2 - The vector V defines a gauge choice. We write V = c\p\ + c 2 p 2 , 
with C\jc 2 ~ We verify that all results are gauge invariant by the 

explicit check that they do not depend on the arbitrary coefficients c\ and 
c 2 - 



3.1 qq — > QQ 

The matrix element for the subprocess qq — ► is very simple. It is given 
exactly by the expression 

(3.6) 

where u and u are the light quark spinors. The equations of motion for the 
initial state quarks has been used to obtain the explicitly gauge invariant ma- 
trix element in eq. ( |3.6| ). The desired Taylor series expansion is obtained by 
using the relation M 2 = 4(m 2 + q 2 ) to expand all factors with M. Converting 
to spherical tensors [O, we find 



M = -^[v(p2h,T a u( Pl )]L^ 



q 2 = -V3[qq}° and a ■ q = -V^Hg, (3.7) 

where the notation [• • -] J M denotes a coupling to angular momentum J and 
helicity M of the spherical tensors inside the square brackets^. The remaining 
Euclidean vectors are converted to spherical tensors after squaring the matrix 
element. 



M s = 2g 2 U c (A j + -e 1 -e 2 z j ) [^a j T c rj- „, „ , - — q'^q^T^). 



3.2 gg -> QQ 

The s-channel gluon exchange diagram can easily be reduced to the form 

h • C2 ^') ~ MW+ 2m) 

(3.8) 

Here e« is the polarisation vector for the initial gluon of momentum pi, and 
T c is a colour generator. For convenience we have used the notation 

Ai = — (ei • Ue 2 ■ pi - e 2 • L»ei • p 2 ) ■ (3.9) 



3 The factor of — \/3 in the conversion of dot products was used earlier in eq. (|2.6|) 
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The t and u channel matrix elements require a little more work. The 
colour factors can be reduced using the identity 

T a T b = U ab + \d abc T c + l -f abc T c . (3.10) 
o 2 2 

The matrix elements can be written in the form 

t u 

M t = 7: and M u = - . (3.11) 

2p ■ pi 2p ■ p 2 

Then the factors of 1/p-pi and \ jp-p2 permit a binomial expansion in powers 
of q. The resulting series in q has coefficients which are simply related to 
t + u and t — u. It is easy to check that 

t + u = - (^p) e Av pM^ &Trj) 

-2M(Ajz m - A m zj + B jm ) (q m ^^r V ) 

+ ( yw + 2m ) S *" B «* (q m q n q p ?° j T V ), (3-12) 
t-u = -2M 2 (A J + ±e l -e 2 z j )(^aJTr ] ) 

Here T stands either for the identity or a generator in the colour SU(3) space, 
depending on which part of the colour structure of eq. ( ^.10| ) we consider. 
We have introduced the additional notation 

Bij = ei ■ Lj62 ■ Lj + 62 • LiE\ ■ Lj. (3.13) 

Note that the binomial expansion is not the desired Taylor series expansion 
in q, since both t + u and t — u involve M, which in turn depends on q. 
However, it is an useful intermediate step, since it allows us to organise the 
terms neatly. 

The full matrix element can be written as 

M = \g 2 5 ab S +\g 2 d abc D c + % -g 2 f abc F c . (3.14) 
o 2 2 

The colour amplitudes S and D involve only Ai t + M-u, whereas F involves 
M s as well as M t — M u . 
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In order to write down our results, we find it convenient to introduce the 
notation 

A = j^e X af,uPiP2 e i4 and <% = A**] + A jZi - Bij + ei • e 2 ZiZj. 

(3-15) 

In order to identify all terms to order v 9 we need the colour amplitude S to 
order g 5 — 

S = - (S«p) A + j_ Sjm {q ^ a j v ) - (mn^ A z m z n (q™q n ? V ) 



M 
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i |p 21 ) Az m z n z p z r (q m q n q p q r ^v) 
M 



S jm z n z p - M + 2 m bjmS np \ z r z s (q m q n q p q r q s g a 3 rj) 

(3-16) 

To all orders, even powers of q come with the tensor A and odd powers with 
S. The amplitude D differs only through having colour octet matrix elements 
in place of the colour singlet ones shown above. For the colour amplitude F 
we need the expansion 

F c = ~ (^?) Az m (q m eT c V ) + ^S om z n (O^V'T^) (3.17) 

In this amplitude odd powers of q come with the tensor A and even powers 
with S. In all three colour amplitudes, the terms in A are spin singlet and 
those in S are spin triplet. 

The decomposition into spherical tensors can be performed partially at 
this stage by using the identities 



Zm,q 



\q\\ and z m z n q m q n = ^-[qq)l - ^[qq)° . (3.18) 



The Euclidean indices on S are most conveniently converted after squaring 
the matrix element. 

A recurrence relation for the i-th term, ti, in either of eqs. ( |3.16|) or ( 3.17|) 
is easy to write. We find that 

U = j^z a z b q a q b ti-2- (3.19) 
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This holds for all i > 3 in the S and D amplitudes and % > 4 for the F 
amplitude. Also, for the F amplitude this holds for i = 3. The % = 4 term 
in F is (2/M)z n q n times the i = 4 term in D. The required Taylor series 
expansion is then obtained by expanding all factors containing M. This 
procedure is completely systematic and may be performed, for example, by 
a Mathematica program. 



3.3 ig -► QQ and 77 -> QQ 



The matrix elements for the two processes jp — ► Q<5 an d 77 
closely related to the amplitudes. It is easy to check that 



M 



19 



geD, 



and 



M 



77 



<5<5 are 



(3.20) 



where D and S are the colour amplitudes given in eq. (|3.14|) , and e is the 
charge of the heavy quark. 



4 The Cross Sections 

4.1 qq -> 

The squared matrix element for this process is easy to write down. After 
summing over initial state helicities, the amplitude square can be expressed 
in terms of matrix elements, over heavy-quark spinors, of products of cr and 
q. At this stage a perturbative projector (eq. |2.8|) is introduced between the 
spinor bilinears in order to project on to QQ final states. The normalisation 
of these states involve the energies of the quarks (see eq. |2.9|) , and can be 
expanded in q 2 , as shown in eq. fl2.10|) . The computation is complete once 



this expansion is performed and the extra factors of q arising from this ap- 
propriately absorbed into the matrix elements. The non-perturbative matrix 
elements needed for the cross sections of various charmonium states are listed 
in Table |1|. 
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H 


d 


Matrix Elements 


Vc 


7 


Os ( 3 S°) 




9 


Os ( 3 S°) 




7 


Os ( 3 S°) 




9 


V s ( 3 S^Sf),O s ( 3 P?) 


Xj 


5 


Os ( 3 S°) 




7 


%( 3 5?, 3 5?) 




9 


8 ( 3 5 1 2 ),P 8 ( 3 5 1 O , 3 ^ 4 ),08 ( 3 ^) 



Table 1: The matrix elements from the qq process contributing to the cross 
section for all charmonium states H at order v d . We use h c as shorthand for 
the l Pi meson. 
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Finally, we list the parton level cross sections- 



cr, 



qq 



* h 

qq 



J/i> 



qq 



Xj 
qq 



R 2 

54m 4 
54m 4 



5{s 



4m 2 )C^( 3 



S(s — 4m 2 



54m 4 
1 

+ — 

mr 



mUV3' 8 



— l—V. 
m 2 1 -y/ - 

5(s — 4m 2 



5?) 

q0 3nz- 
°1> °1, 



1 

— I 

4 



-a J/v, ( 3 pf; 



v^m 1 



12 



12 



(4.1) 

See the appendix for details of the angular momentum coupling scheme used 
in this paper. 

Note that in any application to hadronic collisions, the parton level centre 
of mass energy will be s = X1X2S, where s is the CM energy of the hadrons 
and x\ and X2 are the momentum fractions of the two partons. The con- 
tribution of this sub-process to the hadronic cross section is then obtained 
by convoluting the above cross sections with the appropriate parton density 
functions. 



4.2 gg QQ 

The squared matrix element for the gg process is more complicated, but the 
extraction of the cross section follows exactly the same steps as for the qq 
process. Denoting the average over initial states of the product S* by S ■ S*, 
we find 

S-S* = ^S jm S* /m , = £ [aq}l[ak't x + ^Wq]>' q X- (4-2) 

pol A=±2 Z 

Although a and q are self-adjoint, we have retained the more cumbersome 
notation in order to clarify the coupling of the different angular momenta. 
Also, A - A* = 1/8 and A ■ S* = 0. Consequently, even and odd terms in the 
three colour amplitudes of eq. ( p.l4p do not interfere with each other. 
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H 


d 


Colour amplitude 






S 


D 


F 


Vc 


3 










5 










7 






08 (^l) 


h c 


5 










7 










9 


0i (^o°) 


O 8 (^ 2 ),O 8 ( 3 P]), 


08 {'PD 












J/rI> 


7 




o 8 (^ °),o 8 ( 3 p}) 






9 






08 ( 3 PjO 



Table 2: The matrix elements contributing to the cross section for some 
charmonium states H at order v d . We use h c as shorthand for the 1 P\ meson. 
J = 0, 2 and J' = 1, 2 
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H 


d 


Colour amplitude 


S 


D 


F 


Xo 


5 
7 
9 


V^Pl'Pl) 
o 1 ( 3 P 2 1 ),0 1 ( 3 P 3 ), 

v^pI'p!) 


08(^0°), 


08 CA 1 ), 08 ( 3 ^), 
08 ( 3 ^?) 


Xi 


9 


Oi ( 3 ^ 2 ) 


08 M), 
08 ( S Pi 2 ) 


08 m), 08 ra, 

08 ( 3 ^?, 2 ) 


X2 


5 
7 
9 


01 ("Pi) 

Pi( 3 p 2 Vp|) 

1 (^o 1 ),^ ( 3 P|), 


08 M), 
08 ( 3 4, 2 ) 


08 OP! 1 ), 08 ( 3 ^), 
( 3 ^?,2, 3 ) 



Table 3: The matrix elements contributing to the cross section for some 
charmonium states H at order v d . 
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In Tables ^| and |3| we list all the matrix elements which appear in charmo- 
nium cross sections to order v 9 . The final parton level cross sections are listed 
in the following subsections. The coefficients in the linear combinations of 
matrix elements appearing there depend on the angular momentum coupling 
scheme. Our conventions are set out in the appendix. The contributions of 
these processes to hadronic cross sections are obtained by convoluting the 
sub-process cross sections with appropriate parton densities. 

4.3 Direct J/i/j cross section 

The direct J/ip subprocess cross section is 



a 



J/V, 

99 ' 



s) 



(4.3) 



where 



7r 3 a 2 
4m 2 



5(s — Am 2 



(4.4) 



and ® J J^(n) denotes combinations of non-perturbative matrix elements from 
the colour amplitude a (= S, D or F) at order v n . Using the notation 
explained in the appendix, these can be written as 



(7) 



9^(9) 



e J F ^ 



(9) 



1 



2m 2 
1 



o 



J/i>n 



S° ) + 



2m 4 



v/3 



rrr 



ZO^fP^ + ^fP}) 

_'p J /</73pl 3p3^ 

15m 6 I 4 8 



5 pi 3p3 
r 0' r . 



+2V t 



2m 6 



-a 



J/ip(3pl 3p3 



-o 



J/i>(3p2 



(4.5) 



Previous computations || have considered the expansion only to order 
Heavy-quark spin symmetry gives the relation 0^( 3 Pi) 



upto corrections of order v 2 . Then at order v 7 accuracy this can be used for 
a further reduction of 9p(7). Since we consider the expansion to order t> 9 , 
we cannot use this relation. 
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4.4 xo cross section 

The xo subprocess cross section is 



a 



Xol 
99 y 



^e?(5) + -^e?(7) 

^°(9) + ^e«(9) + A e xo (9) 



(4.6) 



The results to order v 5 are known from previous computations. This work 
identifies all the combinations of matrix elements at orders v 7 and v 9 . Note 
that the order v 7 term involves an off-diagonal operator (i.e., a V term), and 
can be found only after expanding the colour amplitude S to third order in 

Using the notation explained in the appendix, the combinations of non- 
perturbative matrix elements appearing in eq. (|4.6[) are 



e*°(5) 
ef(9) 

0g>(9) 
Qf (9) 



2m 4 

_J^_T>X0(3pl 3p3\ 

1 { °' 0/ 

2 1 

5m 4 1 1 2j + 8m 8 

2m 2 Us 1 ^ oj + 2m 4 



245 
~9~ 



Of( 3 P 3 ) 



4 



149^ 



6m 4 



18m 6 



Of{ z Sl)+hOf{*D\) 



(4.7) 



4.5 Xi cross section 



—ef(9) + —e^(9) + — ef(9) 



(4.8) 



where the combinations of non-perturbative matrix elements are, in the no- 
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tation explained in the appendix, 
1 



&f(9) 

©2(9) 
Of (9) 



2m 4 
1 



W? (3p l) + -Of(3p 2 l) 



08 X1 (^O°) + 



2m 2 8 v uy 2m 4 L 

1 ,a xi ( 1 p 1 1 )+ 1 



Qm 4 



3m 6 



^f( 3 Si 2 ) 



(4.9) 



+ ^8 Xl ( 3 ^)-^8( 3 ^ 2 2 ) 



The Xi is produced first at order v 9 . The large branching ratio for the decay 
Xi — > J/ip makes this a phenomenologically important term, and is the main 
motivation for this work. 



4.6 X2 cross section 



^0F(9) + ^(9) + A©x 2(9) 



where the combinations of non-perturbative matrix elements are 
2 



9* 2 (5) 
Of (7) 

Qf(9) 

©F(9) 

ef(9) 



5m 4 
2 



15m 6 



Of(3P 2 4 ) 



2m 4 
1 



Of( 3 P 1 ) + 



1 



75m 8 
1 



(4.10) 



262 Of ( 3 P 2 3 ) + Hi^pf ( 3 P 2 \ 3 P 2 5 ) 



2m 2 v u/ ' 2m 4 
1 r^8 X2 ( 1 A 1 )+ ' 



9 



2V7 



3m 6 



lores!) + 5 -o?( 3 dd 



-\oreDi) + 2 -on*Di) 



(4.11) 



See the appendix for an explanation of the notation. 
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4.7 rj c cross section 

The production cross section for i] c is- 



-^£(3) + -^6£(5) 



(4.12) 



where the notation in the appendix can be used to write the combinations 
of non-perturbative matrix elements as 



em 

eg (7) 

eg (7) 
eg (7) 



i 



2m 2 
1 



v^m 4 
1 



3m 6 
1 

2m 2 
1 

6m 1 



2(9g(^ 2 ) + ^^( 1 S°, 1 ^ 



(4.13) 



ores*) 

Ot{ l Pl) 



Note that the colour amplitude D first enters the cross section at order v 7 . 
Hence, the almost-elastic cross section 7p — > ?7 C starts at order t> 7 . 



4.8 h c cross section 

The production cross section for the l P\ charmonium state is- 



^(5) + A e ^(7) 



H ^e£(9) + A e £( 9 ) + ^e£(9) 



(4.14) 
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where the combinations of non-perturbative matrix elements can be written, 
using the notation in the appendix, as 



efr(5) 
e&(7) 

©5 C (9) 
©d c (9) 

e h F c (9) 



1 



2m 2 
1 



V3rri 



2O^( 1 S 2 ) + 4EVt( 1 S° n , 1 S^ 
V5 



2m? 
1 

3m 6 



'0; ^0/ 



2m 4 



6m 4 



15m 6 8 v 2y 



(4.15) 



Interestingly, only the colour amplitude D enters the cross section upto order 
v 7 . As a consequence, the cross sections for pp — > h c and 7p — > h c require 
the same combinations of matrix elements upto an accuracy of about 10%. 



4.9 75 -> QQ and 77 -> QQ 

The 7$ cross sections for the production of any quarkonium state, eqs. (|4.3|) - 
( |4. 14| ) , can be obtained from those for the gg process by the following pre- 
scription. Replace a 2 s in cp (eq. |4.4|) by aa s , delete the 0£ and Qp terms, 
and replace the colour factor 5/48 for the terms in G>£> by 2. This follows 
from eq. (|3.20|) . 

Similarly, the 77 cross sections are obtained from the prescription — re- 
place a 2 s in ip (eq. [O]) by a 2 , delete the 0£> and Qp terms in eqs. (|4.3|)-( |4TT4] ), 
and replace the colour factor 1/18 for the terms in G5 by 16. 



5 Discussion 

The number of non-perturbative matrix elements appearing in the cross sec- 
tions in Section § is rather large. In hadron-hadron collisions the influence 
of the qq amplitudes is small compared to that of the gg amplitudes, since 
the gluon luminosity is much larger. This still leaves us with the problem of 
determining a large number of matrix elements. 
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The matrix elements appearing through the colour amplitude D can be 
isolated in almost-elastic 7p collisions. At HERA it is now possible to sepa- 
rate the diffractive components of the cross section |T2[. With this separation, 



it becomes possible to measure these matrix elements with greater accuracy 
than was possible in the past [[0| . 



The matrix elements arising from the colour amplitude S can be mea- 
sured in 77 collisions produced in e + e~ colliders, provided the photons are 
point-like. Unfortunately, cross sections for 77 collisions at the LEP are 



likely to be dominated by resolved-photon effects |T/J. On the other hand, 



double-tagged events with low hadronic energy at the LEP correspond to 
the scattering of two highly virtual photons, 7*7*. The cross sections for 
7*7* — > J ft/) (and other quarkonia) can also be computed in NRQCD. It 
is easy to see that these cross sections are obtained from the colour singlet 
part of the amplitudes in eq. ( ft.llj ), the only difference being in the denom- 
inators of the propagators. Hence, all the matrix elements from the colour 
amplitude S which appear in Section |4] will also appear in this case, albeit in 
different linear combinations. Furthermore, the coefficients would depend on 
the virtuality of each photon, leaving us with the possibility of independent 
measurements of each of these matrix elements. The full computation for 
such 7*7* processes will be presented elsewhere. 

The colour amplitude F cannot be separated in any process, since it arises 
only through gg initial states. The corresponding matrix elements may be 
obtained from pp, np, pp collisions, as well as in 77 collisions where the 
resolved photon processes dominate. 

Although quantitative estimates of low transverse momentum quarko- 
nium cross sections involves the prior analysis of many such experiments, it 
is possible to make some crude but interesting estimates. 

The gg — > g c cross section (eq. (4.12| ) involves only the colour amplitude S 
to next-to-leading order, v 5 . The colour amplitude D first enters at order v 7 . 
As a result, the 7p — > r\ c cross section must be much less than the pp — > r\ c 
cross section. In contrast, only the colour amplitude D enters the gg — ► h c 
cross section upto the next-to-leading order, v 7 . Thus the pp — > h c cross 
section can be predicted to better than 10% accuracy if the ■yp — > h c cross 



section is known. If the 1 P\ charmonium state is identified definitively [15 
then this fact could serve as an empirical test of NRQCD factorisation. 

Further qualitative arguments depend on a more detailed analysis of the 
non-perturbative matrix elements. Any such matrix element O, with mass 
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dimension D, may be written in QCD as 



O = A D f (^) , (5.1) 

where A is the QCD scale and rrif are the masses of the quarks of flavour /. 
In the chiral limit, the light quark masses may be taken to vanish, and the 
corresponding arguments of / are then zero. For the matrix elements which 
appear in quarkonium cross sections, the dependence on the heavy quark 
mass is factored into the coefficient functions. After such a factorisation, 
the resulting matrix element for the production of a quarkonium H can be 
written as 

O* = A D f H (v), (5.2) 

where v is the dimensionless velocity which organises the NRQCD expansion. 
In the limit v — > 0, f H — > c H v d , where c H is a constant and d is the power 
counting dimension obtained, for example, by eq. ( |2.5|) . Power counting in 
NRQCD is reliable if and only if c H obtained from different operators (but 
for the same H) are of similiar orders of magnitude. 

Heavy-quark spin symmetry gives rise to further restrictions on these con- 
stants, c H . The operator expectation values from different hadrons may be 
related using this approximate symmetry. Then, for the same operator, the 
values of the constant obtained with two different hadrons with the same ra- 
dial wave-function, must be also of similiar magnitude, after removing certain 
Clebsch-Gordan coefficients. 

Such arguments may be used to study the convergence of the NRQCD 
series for various cross sections. In the approximation of eq. ( p.2[ ), the coef- 
ficient function times the matrix element is a function of two dimensionless 
variables, v and z = A 2 /m 2 . With the assumption that all the constants c H 
are of similiar order, the convergence properties are simple, since z<d 2 . As 
an example, we take the cross section for gg — > r] c (eq. 4.12|) . After removing 



all common factors, including dimensional quantities, we find 

z 2 + 0(z*). (5.3) 



99 



—v 3 + —v 7 + 
36 96 



The first term comes from O v s c {3) and the second from 0^(7) in eq. (|4.13|) . 



The remaining terms are of higher order in z, and converge rapidly. The two 
coefficients in eq. ( |5.3| ) are of similiar order, and it is quite possible that the 
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series in v 2 may be convergent for small v. Nevertheless, it is interesting to 
note that the main correction to the leading term comes, not from the order 
v 5 terms, but from one of the order v 7 terms. 

A similiar argument may be used to assert that the most important term 
in the cross section for gg — > Xj is the 0£>(9) term (eqs. [477] , [O] and |4.11|) , 
since it has the operator with the lowest mass dimension. In that case, 
heavy-quark spin symmetry may be used to show that 

a( Xo ) ■ <t(xi) ■ vfa) ~ 1:3:5, (5.4) 

both in pp and almost-elastic jp collisions. The measured value of this ratio 
in pp and Tip collisions is 0.62±0.18±0.09 [jTE| . Corrections to eq. ( |5.4j ) would 



come from the breaking of the approximate heavy-quark symmetry at higher 
orders in v. Consequently, the same ratios should also be observed for the 
corresponding bottomonium states. Since the coefficients for quarkonia with 
unequal principal quantum numbers cannot be compared, the prediction of 
the direct to total J / ip cross sections is a more detailed dynamical question. 

It is interesting to test such a scaling hypothesis against the matrix ele- 
ments which have already been extracted. Unfortunately, a definitive test 
cannot yet be made, since the values quoted for the same matrix elements 
vary widely |T7], [TJ| . We take the examples of the two sets 

(<?f( 3 Si)> and \(0»{ 1 S Q )) + — 2 (O^Po)), (5-5) 

for J/ip, ip', T(1S) and T(2S). The two L = matrix elements give a 
common coefficient c H = ci, and the L = 1 matrix element gives a differ- 
ent coefficient c-i- For the same H, if the scaling ideas are correct, then 
(ci — Oi)/(ci + C2) should be approximately zero. We find that for the T 
states the fitted numbers [17|| are within one standard deviation of this re- 



sult. The values for ijj' [O are about two standard deviations away, whereas 



the many different fits for the J/ip [0, [18| lie between five and two standard 



deviations from this expectation. Although one is tempted the accept the 
scaling argument, it turns out that c 2 is systematically larger than c\. This 
clearly calls for an inclusion of higher order effects, both in a s and in v. The 
latter, specially, may give rise to unforeseen corrections in view of the pattern 
of coefficient functions seen here. 

We would like to thank V. Ravishankar for discussions. 
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A The Coupling Scheme 



We use spherical tensor techniques for the reduction of the matrix elements. 
Since all these are constructed by multiple couplings of Euclidean 3-vectors, 
we need the spherical tensor components for any vector V — 

V = V g , Vx^-^=(V x + iV y ), V-i = -^=(V x -iV y ). (A.l) 

The labels on the spherical tensor components denote helicity. Couplings 
of spherical tensors require the usual Clebsch-Gordan coefficients fll!] . We 
denote coupled tensors by the notation [a, b] J M , where J is the rank of the 
coupled tensor and M is the helicity 

Next we consider the complete specification of the state /Q|0), where /Q is 
one of the operators in eq. (|2.1| ). When computing the NRQCD cross section 
to high orders in v, we may have to couple a large number, N, of g's, and, 
possibly, a Pauli matrix. It is possible to write down the Clebsch-Gordan 
series for a reduction to states of given J, J z , L and S, keeping explicitly the 
permutation symmetry of the g's. However, the coupling constants are not 
listed in the literature, except for N = 3. We choose instead to first ignore 
the permutation symmetry, treat all the g's as distinguishable operators, 
and use the usual methods of angular momentum recouplings. When this is 
completed, we explicitly symmetrise the results to get the desired expressions. 

The total number of commuting operators required to specify the state 
in the direct product basis is 2N + 2S. In the coupled basis we have N labels 
from each derivative operator, 2 from the total angular momentum J and 
the helicity J z . If S ^ then two more labels, L and S have to be specified. 
All these total to + 25* + 2. Consequently another N — 2 labels have to be 
given for a complete specification of the state. It is sufficient for this purpose 
to fix an order of coupling the individual derivatives and the Pauli matrix 
and specify the tensor rank at each coupling |TU . 



As an illustration we take the case S = 1 and N = A. We specify a 
coupling order either by the bracketed expression on the left or the binary 
tree on the right — 



a 



[qq] n [qq] 



32 33 



.1 



(A.2) 



Each node in the tree denotes a tensor whose angular momentum is to be 
specified. The root node (at the bottom) is the total angular momentum 
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J. The helicity, J z , is also specified at this node. All the leaf nodes (at 
the top) correspond to the individual Euclidean vectors being coupled. A 
circled node represents a Pauli matrix. A simple counting shows that all the 
required labels for the state are specified by this means. 

In the Taylor expansion of the appropriate matrix elements (section ||) 
there is a natural pairing of derivatives. The factors of z m z n in the recurrence 
relation eq. ( [3.19| ) naturally pair these operators. In addition, 6 mn from q 2 
factors or the tensor S np in the fourth term of eq. ( |3.16|) also pair derivatives. 
The tensor Sj m pairs a Pauli matrix with a derivative. In addition, there is an 
unpaired derivative in the qq amplitude (eq. |3.6| ) and the colour amplitude F 
(eq. |3.17|) . A natural set of groupings is then given by the following scheme — 

1. The spin singlet terms in the colour amplitudes S and D have even 
number of derivatives. These are grouped as 



[qq] ■ ■ ■ [qq] 



(A.3) 



The spin singlet terms in the colour amplitude F have an odd number 
of derivatives. Exactly one is unpaired. The rest are first coupled, and 
this unpaired derivative is coupled at the end — 



[qq] ■ ■ ■ [qq] 



(A.4) 



The spin triplet terms in the colour amplitude F or the qq amplitude 
have an even number of derivatives. Exactly one is unpaired, and 
one paired with the Pauli matrix. The [qq] pairs are first coupled; 
the unpaired derivative is coupled to this, and finally the [erg] pair is 
recoupled to the result using a 6 — J symbol. A further recoupling with 
another 6 — J symbol then joins the two uncoupled derivatives — 



a 



[qq] ■ ■ ■ [qq] 



a 



[qq] 



[qq] ■ ■ ■ [qq] 



(A.5) 



4. 



The spin triplet terms in the colour amplitudes S and D have an odd 
number of derivatives. One is paired with a Pauli matrix. The rest of 
the derivatives are first coupled as before, and the remaining derivative 
is then recoupled to this. The coupling scheme is 



a 



[qq] ■ • • [qq] 



(A.6) 
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In this work N < 5, and consequently, a choice of the order of couplings of 
more than two pairs [qq] is not required. At higher orders further choices 
have to be made. Except at low orders (N < 2), the choice of scheme is not 
unique. Different coupling schemes give rise to different sets of basis states 
for /Q|0). Unitary transformations can always be found to relate these sets to 
each other. Although the linear combinations of matrix elements appearing 
in the cross sections then look different, they can be transformed into each 
other. 

A possible phase ambiguity might remain in the definition of the states 
because the coupling orders [a, 5p and [b, ap differ by the phase (— . 
Enumeration of all the cases arising in this problem shows that there is no 
such ambiguity — 

1. Since [qqY has only j = and 2, the phase is 1, and no ordering 
ambiguity exists. 

2. For [[w]a[<?<?]o]a' if A = 0, then we must have k + j + 1 even. For integer 
k and j the phase is then 1, and no ambiguity exists. When A = ±2, 
the phase is 1 if j = 0. If j ^ 0, then the phase is not identically 
1. However, in this case we always have a sum over helicities with 
multiplicative factors which are independent of A. This projects out 
the terms with k + j + I even, and removes the ambiguity. 

3. For [[g]o[[<?g] • • • [99]]a1a> the ordering is immaterial if A = 0, since 1+1+ j 
is then even. When A = ±2, the sum over A removes the ambiguity. 

4. For [ergp, since j = or 2, the phase is 1. 

5. For the final L — S coupling, we always choose the Pauli matrix as the 
first factor. 

The binary tree is thus a complete and unambiguous specification of the 
coupling scheme for our purposes. 

After the recouplings are completed, we recognise that the tensors are 
symmetric under interchange of all g's. A symmetric tensor with N 3-d 
Euclidean indices reduces under the rotation group. The allowed angular 
momenta are L = N, N — 2, etc, down to or 1, depending on whether N is 
even or odd. Simply counting the number of components in the tensor shows 
that each allowed value of L appears only once. As a result, it is sufficient 
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to label the symmetric states by L and iV instead of all the intermediate 
angular momenta. 

The unitary transformation between the scheme specified by all the in- 
termediate angular momenta and the symmetric tensors can be carried out 
explicitly. We find that — 



M 



aiqiqq] ] 1 




a 



a 



a 



a 



q[[qq]°[qq] oVI 
q[[qq] 2 [qq] 2 ]° 

q[[qq)°[qqF 2 



q[[qq] 2 [qq] 



2l2' 



■ 1" 


J / 




M \ 


- 1" 


J j 




M \ 


■ 1" 


J I 




M \ 


■ 1" 


' - ( 




M \ 



a \q\qq} 2 } 1 
[[?9] 2 [??] 2 ]S 

i Ff 



j 

M 



3p3 

3/ J 

1 c4 

On 




(A.7) 



V3Vi5y 



3 p5 



We specify the state /Q|0) by a small extension of the spectroscopic no- 
tation. In most cases we will write the state, or equivalently, the operator, 
as 2S+1 Lj , where gives the number of derivatives. This allows us to use 
the new index A^ for power counting (see eq. |2~H| ). 

The matrix element (K.iU(H)JCj) is clearly zero unless common quantum 
numbers of the two states /Cj|0) and /Cj-|0) agree. Thus L, S, J and J z 
must be the same. Subject to these constraints, no quantum number reason 
prevents a non-zero matrix element for bilinears with unequal N. Hence we 
retain such off-diagonal operators in our computations. 
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